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i\  „ rN'i'Romrcn  on 

Formulas  for  the  computation  of  the  electric  surface  current  and  far 
.scattered  field  of  a perfectly  conducting  body  of  revolution  are  derived 
for  arbitrary  plane  wave  excitation.  Computer  program  subroutines  which 
implement  these  formulas  in  the  resonance  region  will  appear  in  a forth- 
coming report.  Computations  show  that  both  the  H-field  solution  and  the 
E-field  solution  deteriorate  near  internal  resonances  of  the  conducting 
surface,  but  that  the  combined  field  solution  does  not. 

The  field  solutions  are  obtained  by  applying  the  method  of  moments 
to  the  H-field  integral  equation,  the  E-field  integral  equation,  and  the 
combined  field  integral  equation  for  a perfectly  conducting  body  of  revo- 
lution. Although  the  computer  program  subroutines  are  written  explicitly 
for  a perfectly  conducting  body  of  revolution,  they  are  directly  applica- 
ble to  the  more  difficult  problem  of  plane  wave  scattering  by  dielectric 
bodies  considered  by  Wu  [1]. 

Our  H-field  solution  is  similar  to  that  of  Uslenghi  [2],  generalized 
to  oblique  incidence  and  with  expansion  and  testing  functions  equal  to 
four  impulse  approximations  to  triangle  functions  divided  by  the  cylin- 
drical coordinate  radius.  We  use  Gaussian  quadrature  instead  of  Simpson's 
rule  for  integration.  Our  treatment  of  coincident  impulses  is  simpler 
than  Uslenghi 's.  The  impulses  are  combined  in  groups  of  four,  not  to 
make  the  H-field  solution  more  efficient,  but  to  make  it  compatible  with 
the  E-field  solution.  Actually,  for  low  order  solutions  where  much  more 
effort  is  required  to  obtain  the  matrix  elements  than  to  solve  the  system 
of  linear  equations,  it  is  wasteful  to  combine  the  impulses  if  all  one 
wants  is  the  H-field  solution  alone. 

Our  E-field  solution  is  in  some  respects  a simplification  and 
in  other  respects  a refinement  of  an  earlier  E-field  solution  [3). 


[1]  T.  K.  Wu,  "Electromagnetic  Scattering,  from  Arbitrarily-Shaped  Lossy 
Dielectric:  bodies,"  I’li.I).  Hies  is.  University  of  Mississippi,  May  1976. 

[2]  P.L.E.  Uslenghi,  "Computation  of  Surface  Currents  on  Bodies  of  Revolu- 
tion," Alta  Kretpienza , vol  . 39,  No.  B,  1970,  pp . 1-12. 

| 3 1 .1 . it.  Maut  z and  K.  U.  Harringt  on,  "Radiat  ion  and  Scat  tering  f rom  Bodies 
of  Revolution,"  App  I . S e i . _Res  . , vol.  20,  June  1969,  pp  . 40.S-4  1S. 
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The  interaction  between  impulse  portions  of  the  expansion  and  testing 
functions  is  calculated  in  the  same  way  as  in  the  H-field  solution. 
Computationally,  this  new  E-field  solution  is  roughly  three  times  faster 
than  that  of  [3]  with  comparable  accuracy. 

Our  combined  field  formulation  is  that  proposed  by  Oshiro  et  al. 
[4,5].  It  is  obtained  by  the  method  of  moments  applied  to  a weighted 
average  of  the  H and  E field  Integral  equations  for  a perfectly  conduct- 
ing body  of  revolution.  The  matrix  operator  for  the  combined  field 
solution  is  a linear  combination  of  the  matrix  operators  for  the  H and 
E field  solutions.  The  excitation  vector  for  the  combined  field  solu- 
tion is  the  same  linear  combination  of  the  excitation  vectors  for  the 
H and  E field  solutions. 


II . STATEMENT  OE  THE  P RDRLF.M 

We  seek  the  electric  surface  current  and  the  far  scattered  field 
of  the  perfectly  conducting  body  of  revolution  of  Eig,  1 excited  by  an 
incident  plane  wave.  In  Fig.  1,  p,i,z  are  cylindrical  coordinates, 
and  t,4>  form  an  orthogonal  curvilinear  coordinate  system  on  the  surface  S 

of  the  body  of  revolution.  Also,  u.  and  u are  orthogonal  unit  vectors  in 

"*  t.  '“\j> 

the  t and  <f,  directions,  respectively.  The  coordinate  origin  is  on  the 

axis  of  the  body  of  revolution  hut  not  necessarily  at  the  lower  pole  a^ 

in  Fig.  1,  Figure  2 defines  the  propagation  vector  k of  the  incident 

plane  wave,  the  transmitter  coordinate  6 , the  coordinates  6 , d>  (re- 

t r r 

ceiver  coordinates)  at  which  the  far  scattered  field  is  observed,  and 

the  propagation  vector  k of  a hypothetical  measurement  plane  wave  which 

~r 

travels  from  the  receiver  location  (fl  , <f>  ) toward  the  origin.  Note  that 
the  1 coordinate  of  the  transmitter  is  zero  such  that  k is  in  the  xz 

v~t 

t t r r 

plane.  Tn  Fig.  2,  uA  , u . , u.,  and  u,  are  unit  vectors  in  the  0 , y,  0 , 

’*•{)  A'kH  -^f)  £ y 

and  (Ji  directions  respectively. 


We  consider  separately  a 0 polarized  incident  plane  wave  defined 


[4]  F.K.  Oshiro,  K.M.  MItzner,  and  S.S.  hocus  et  al.,  "Calculation  of  Radar 
Cross  Section,"  Air  Force  Avionics  Laboratory  Tech.  Kept  . AFAI.-TR-70-2  I , 
Port  II,  April  ]«70. 

[5]  A.  !.  Pog.gio  and  F.K.  Miller,  "Integral  Kquat  ion  Solutions  of  Three- 
Dimensional  Scattering  Problems,"  Chap.  4 of  Computer  Techniques  for 
Electromagnetics,  edited  by  R.  Mittra,  I’ergainon  Press,  1073. 


Body  of  revolution  and  coordinate  system. 


'Mane  wave  scattering  by  a conducting  body  of 
revol ut i on . 


by 
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and  a <j>  polarized  incident  plane  wave  defined  by 


(1) 


f-  -jk  * r 

E = u j kp  e 

*■<!> 


--jk  • r 

,.i  ,.  t.  •'  »t 

H = (k  x u )e 

V.  .V  £ V*-  (J} 


(2) 


where  E1  and  H1  denote  incident  electric  and  magnetic  fields  respec- 
tively, r is  the  radius  vector  from  the  origin,  k is  the  propagation 
constant,  and  n is  the  intrinsic:  impedance.  Either  plane  wave  gives 

rise  to  t and  cf>  directed  electric  surface  currents  on  S and  0 and  6 

r r 

directed  far  scattered  fields. 


HI  . H- FIELD  SOLUTION 

The  H-field  solution  is  obtained  hv  applving  the  method  of 
moments  to  the  Il-field  integral  equation.  The  H-field  integral  equa- 
tion is  derived  by  setting  the  component  tangential  to  S of  the  total 
magnetic  field  equal  to  zero  just  inside  S. 

'Hie  boundary  condition  that  the  total  tangential  magnetic  field 
is  zero  just  inside  S is  written  as 

-n  * HS  = n * li1  just  inside  S (3) 

where  n is  the  unit  outward  normal  vector  to  S,  llK  is  the  magnetic- 
field  due  to  the  electric  surface  current  on  S and  li'  is  the1  incident 
magnetic  field  given  by  either  (1)  or  (2). 

To  obtain  an  expression  for  n if’,  we  note  first  that  from 
page  98  of  [6j 


(6)  K.  i'.  Harrington,  Time-Harmonic  Electromagnetic  Fields,  Mi  < lraw-ll  i 1 1 
Book  Co . , 1961. 
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to 


lim  n x V x 


S 


(4) 


n x H 
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-ikin' | 

.T(r')  2 ds* 

4 it  | r-r ' | 


where  6 is  the  distance  be.tween  the  field  point  r and  the  inside  face 
of  the  surface  S,  r'  is  a running  source  point  on  S,  and  J(r!)  is  the 
electric  surface  current  on  S.  Next,  we  view  J(r')  in  (4)  as  the 
current  which  resides  on  AS  plus  the.  current  on  S minus  AS  where  AS 
is  that  portion  of  S inside  a sphere  of  radius  e centered  at  the  point 
on  S nearest  r.  Let  c be  so  small  that  AS  is  essentially  flat  and  that 
the  electric  surface  current  on  it  is  constant.  If  d is  appreciably  less 

g 

than  r , then,  from  Ampere's  law,  the  contribution  to  n * H from  the  cur- 
rent on  AS  is  -T/2  where  J is  the  value  of  the  current  on  AS.  Moreover, 
this  contribution  to  nx  IIs  comes  exclusively  from  a small  portion  of  AS 
in  the  immediate  vicinity  of  r.  The  current  on  that  portion  of  AS  for 
which  the  distance  to  the  field  point  r is  appreciably  greater  than  A does 
not  contribute  to  n*Hs.  We  now  let  6 -►  0 in  which  case  r becomes  a point 
of  S,  the  contribution  -.1/2  comes  from  the  value  of  J at  the  single  point 
r,  and  the  current  on  any  portion  of  AS  which  does  not  contain  r e.ontri- 

g 

butes  nothing  to  n * 11  . Hence, 

s I(I}  ff  e-lklrT-'l 

n * If  - - -----  + II  n x Vv  tJ(r')  — ]ds’  (5) 

' " 2 jsj  - " 4,  | r-r*  | 

wlu're  r is  exactly  on  S and  where  the  improper  integral  In  (5)  is  con- 
vergent . 

In  view  of  (5)  and  the  vector  idem  i. ty 

• I ’.(>"•) 

t i)  becomes 


j k | r-r 1 j If  j k | r-r ' ! 
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(- 


. — i k [ r-r' 

— )e  v. 


n x 


[(r-r')  x ,T(r ' ) ] ds  ' 


x H (r) 


r-r 


for  r on  S . Henceforth,  we  assume  that  the  outward  normal  vector  n is 
Riven  bv  n = u x u . if  u is  chosen  such  that  n = - ti * u , then  our 

■ v,  -t  ^ ~ t 

evaluations  of  the  terms  proportional  to  n in  (7)  will  have  the  wrong 
signs  throughout  the  remainder  of  this  report. 

If 


T(r')  = u\  Tt(t\<fr,>  + u’  (8) 

^ +*  f;  ^9 

where  u'  and  u'  are  unit  vectors  in  the  t'  and  * directions  respec- 
tively,  then,  as  shown  in  Appendix  A,  (7)  becomes 


i-Vjj-  + | p'dt'  j d.(!  '(kTC  (t  ',  [ ( (p  *-p)  cos  v'  -(z  ' -z)sin  v') 


0 


. 3 r 


cos  !>'-2p  cos  v's  in"  (~)  ] -t-  — j p'dt'  j d 1 '(1TS  (t ' , |>  '+<|>)  (z  ' ~z)  sin  1 ' 


+ J i ' (t , i)  f k3  r 

+ 'U}~~2 4-r  J p llt 


I ■! ' fl.T  ^ (t  ',*'  + !)  (p  's  in  v cos  v'-psin  v'cos  v 


( z ' ) sin  v sin  v ' ) sin  ; ' + -- - ! p 1 <1 1 ' 

‘i  ■■ 


! 

<*’!  i 


d-f'dff  (t  ' '+,•) 


> t ' . 


[((,  '-p)cos  v - (z'-z)sin  v)  cos  ;> ' f 2,  Vos  v sin"  (-“  ■)  1}  n - II  (9) 


who  re 


I + j k K - 1 k K 
. - - .■ 

k K 


(nn 


n 


(11) 


R = V(p--i  i')^  + (z-*z')*"  + 4np'  s i n ‘ (•'-■-) 


In  (9),  both  n and  Hl  art-  to  be  evaluated  at  (l,'!0  on  S and  v Is  the 
angle  between  u and  the  v,  axis.  v is  positive  11  u points  away 
from  the  z axis  and  v is  negative  if  u points  toward  the  z axis.  The 
variables  p ' , z',  and  are  respectively  p,z,  and  v evaluated  at  t ’ . 
If  the  electric  surface  current  is  bounded  and  if  S curves  smoothly, 
then  all  the  iterated  integrals  in  (9)  converge  because  the  Integrands 
are  at  least,  as  well  behaved  as  R ^ . 

According  to  the  method  of  moments,  we  let 


.T(r)  = V (l1.  Jt.(t,^)  + I*.  J ^ . ( t , <j>) ) 
- ~ e . ni  v-n  t ni  - n -t  ^ 

n , 1 J 


02) 


t (j) 

■here  J . (t,<f)  and  J . (t,<f>)  are  expansion  functions  defined 

■~n  t A-n  i 


u-tri(t) 


e:in«J) 


(13) 


In  J. 

e 


(14) 


Expansion  functions  whose  t and  1 components  are  proportional  to 

e-in*  are  especially  suitable  because  they  make  the  t and  f components 

i n (ti 

of  the  left-hand  side  of  (9)  proportional  to  e ' . The  coefficients 

T . and  T^.  are  determined  bv  solving  the  matrix  equation  which  results 
n.l  n | 

when  (12)  is  substituted  via  (8)  into  (9)  and  the  inner  product  of  (9) 

t 1 

with  testing  functions  W ^ ( t , <Jj ) and  W_‘  (t,f<)  defined  bv 


WtJ  = u f,(t)  .*  (15) 

-mi  - t i 

W ' . = u f (t ) e ^ ' (lb) 

mi  ! 
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is  taken.  For  an  inner  product  equal  to  the  dot  product  integrated 
over  S,  the  matrix  equation  decomposes  into 


’ y“ 

n 

ytp 

n 

" r 

n 
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n 

Y^ 

i* 

fib-i 

n 

n 

n 

- n I 

, n=0 , +1,  +2. . 


(17) 


~ (f)  £ 

where  I and  I are  column  vectors  whose  i th  elements  a~e  I and 
, n n v . , ~ n i 

(f*  r 3~d) 

1 . respectively.  Also,  1 and  Tr  are  column  vectors  whose  ith  ele- 
ni  n n 

ments  are  given  hy 


f . = dtpf.(t) 
ni  i 


and 


ni 


dt  pf.(t) 


2-'T 
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2 IT 
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• H1 


d<J>  (u  k n)  • H e 


d<f>  (u  x n)  • H1  e 3 n ' 

"<t>  ■' 


(18) 


(19) 


respectively.  Finally,  Ytt;,  Ylt:,  Yt3,  and  are  square  matrices 

nan  n 

whose  ij th  elements  are  given  by 


(YtC)  = tt  f dt  p f.(t)f.(t)  !-k3  I dtpf.(t)  dt’p’f.  (t 1 ) [(Cp'-p)cos  v' 
^ V 1 J 1 .1  J ^ J .1 


(z*-z)  sin  v')  “ r,  p cos  v' 


(20) 


(Y'Pt)..  = ik3 

n ij 


dt p f (t) 


d t ' p ' f . ( t 1 ) (p'sin  v cos  v'  -P  qln  v'  cos  v 


- (z'-z)  sin  v sin  v')C, 


(21) 


(Y^  f ) . . = j k 3 | dt 


n i 


> f.(t)  f dt'p'f.(t  ’)(z’~z)(, 

! J 1 2 


(22) 


fi 


' ) . . = n j dtp  f , ( f ) f . ( t ) + k ^ 


dtp  ft(t)  j dt’p'f.d ’)[((p*-p)(,os  v 


(z'-z)  s in  v)G0  + Op’  cos  v] 


(23) 


where 


r.  = 2 


d(j)’G  sin"  (4> ' /2)  cos  (nJ>') 


(24) 


C.  = 


d^'G  cos  ip " cos  (mb') 


(25) 


0^  = j dl'G  sin  <p v sin  (ivt1') 
0 


(26) 


Wo.  define  pf  (t)  to  be  a four  impulse  approximation  to  a triangle 
function  in  the  following  manner.  Letting  t = (p,z)  denote  that  p and 
z are  cylindrical  coordinates  of  the  point  t,  we  define  an  odd  number 
greater  than  or  equal  to  5 of  consecutive  points  t = (p,z)  = t^  = 

(p^,  Zj ) , i=l,2,...P  on  the  generating  curve  of  the  body  of  revolution 
such  that  (pj,  z^)  and  (pp,  Zp)  are  the  poles.  If  the  body  of  revolu- 
tion has  no  poles  because  the  generating  curve  closes  upon  itself  as 
with  a torus,  then  three  points  must  be  overlapped  such  that 


( h P _ 


3+i! 


P-3+i 


.) 


V V- 


1,2,3  . 


Preferably,  the  points  t.  should  be  such  that  n 


u,  * u where  u 
-(J>  -t  ~t 


points  from  t.  t; 


i + l 


f c the  points  t are  chosen  such  that 
i 


n - - u^ut,  then  all  expressions  which  can  be  traced  hack  to  the 
terms  proport ional  to  n in  (7)  will  have  the  wrong,  signs. 


We  now  approximate  the  generating  curve  by  drawing  straight:  lines 
between  the  points  (p  ,z^),  1=1, 2,... P and  define  points 


4 


, N /’i  + Pl+1  "’I  ‘ :i14ls  ,.„-,N 

t ^ - (P ! , *,>  = ( 2 * 2 } U/) 

on  this  approximate  generating  curve.  The  length,  d of  the  interval 
centered  about  t is  given  by 

ji  - /<*;«  - »i>2  + (:t+1  - v2  «» 

In  terms  of  coefficients  T ^ ^ defined  by 

= kd2i-l 

4i-3  2'd2i-i  + d2i‘> 

T = k(d2i-l  + 2 d2.i)d2i  (29) 

4i-2  d^  ^ + d2. 

= k(d2i+r.  + 2 d2i+l)d2i+l 

4i-J  d„ . . , + d . 

Z lfl  2 1+2 

tJW2  _ 

4i  - 2(d21+1  + d,.+2)  • 

we  construct 

’V1*  ’ i I,  V'.i-4  s<t-tP+2i-2)  °0> 

P=1 

wliere  4(t)  is  the  unit  impulse  function.  The  right-hand  side  of  (30) 
in  the  desired  four  impulse  approximation  to  a triangle  function 
(see  Fig.  3). 

Substitution  of  (30)  into  (20)  - (23)  yields 


10 


Ik 


(y“)  . 

(Yl) 

n ij 

s 

(Ynt}i.1 

4 4 

(Y2)„ 

b 

(Yt<f\. 

= y i’  y 7' 

p+4i-4  q+4i-4 

p = l K q = l M ■’ 

(Y3) 

n i.l 

s 

(Y*V. 

(Y4) 

. n i.l  a 

s 

where  s denotes  the  double  subscript  p+2i-2,  q+2j-2  and 


(Yl) 


*■1  ) 


k((p.~p.)  cos  v.  - (z,-z.)  sin  v . ) C„-kp . G.  cos  v.,  i^j 
1 i .1  .1  i 1 2 l 1 l 


(32) 


l 


— x - kp  . G cos  v. 

,2  i 1 i 

k cl . p . 
l l 


> 1=1 


(Y2)  . . = ik(p.sin  v.cos  v.  - p.sin  v.c.os  v. -(z  . -z  , )sin  v.sin  v.)G„  (33) 
ij  1 i ii  1 i 1 i i 1 3 


<Y3>11  ‘ 1k<2,  - 2i)G3 


(34) 


k((p.-p.)cos  v.  - (z.-z.)  sin  v.)G„+kp .G, cos  v.  , i^i 
I i i j i l 2 ] 1 l 


(Y4)i.i  = 


(35) 


k2diPt 


+ kp.G,  cos  v. 
i 1 l 


i=l 


Here,  is  the  angle  between  the  approximate  generating  curve  at 

(p.,z.)  and  the  z axis.  The  term  ■ — — in  (12)  and  (15)  was  obtained 

1 1 k d.p. 

l < 

by  replacing  one  of  the  coincident  impulse  functions  in  t lie  first  inte- 
gral on  the  right-hand  sides  of  both  (20)  and  (21)  by  an  equivalent  pulse 


over  the  interval  of  length  d . 
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In  (,  J2)  ( 55)  , C^,  G^ , and  G^  are  given  by 
C is  given  by  (10)  with  R of  (11)  evaluated  at  (p 
When  i = j in  (32)-(35),  none  of  the  integrals  G 
cause 


(24)  to  (26)  in  whiclt 


,z,p  ' ,Z  ')  = (p  . ,7  ,p 

,0^,  and  G converge 


R = 2p^sin  (if  ' /2) 

This  lack  of  convergence  is  ascribed  to  use  of  the  impulse  rep re sent a 
tion  (30)  rather  than  to  any  deficiency  in  the  H-field  integral,  equa- 
tion (9).  To  obtain  convergence  for  i = j,  we  replace  the  above  R by 
an  equivalent  distance  Rg  given  by 


R 

e 


= ■/  (d./4)2  -I-  4P2sin2(<f,72) 


(36) 


Expression  (36)  may  be  obtained  by  displacing  the  field  point  a distance 
di/4  perpendicular  to  the  plane  of  the  source  loop.  Now,  d /4  is  the 
equivalent  radius  [7]  of  a flat  strip  whose  width  is  d_.  . Expression  (36) 
can  also  be  obtained  by  averaging  R2  for  field  points  displaced  a dis- 
tance d±/A  in  either  direction  along  the  approximate  generating  curve 
from  the  source  loop. 


An  point  Gaussian  quadrature  formula  is  used  to  calculate  the 
integrals  G^,G^rand  defined  by  (24) -(26)  in  which  G is  given  by  (10) 
and  R by  either  (11)  or  (36).  According  to  this  quadrature  formula. 


f ( ' ) d (J)  ’ = 


~ r 

i l 

k.=l 


\r(f  <\  + >» 


(37) 


where  f(<j»')  is  the  function  being  integrated  and  and  A are  constants 
tabulated  by  Krylov  [8].  In  (37),  the  multiplier  ~ and  argument  ~(x  +1) 


17] 

[8] 


E . A. 

1 966 , 

V.  T. 

A.  II. 


Wolff,  / Jitenna  An  a iy  s i s , dohn  Wiley  and  Sons,  Inc. , New  York 

p . 61  . 

Krylov,  App vox i mate  Galculation  of  Integrals,  translated  by 
Stroud,  Macmillan  Co.,  New  York,  1 962*/  Append ix  A. 


1 3 


A 


instead  ol  just  x are  due  to  t lie  t ransformat  ion  of  Krylov’s  interval 
from  -1  to  1 into  the  interval  from  0 to  n . 

The  — terms  in  (32)  arid  (35)  destroy  some  symmetry  prope r- 

k "d  . n . 
x r 

7T 

ties  of  (31).  If  the  — - — terms  were  absent,  then 

k“d . p . 

l i 


(YI)  . . 
i.l 

' - (Mhi 

(Y2)  . . 
11 

- - <V2)U 

(38) 

(YDii 

" - CV3) .. 

.11 

and  replacement  of  (i,j,p,q)  by  (j,i,q,p)  in  (31)  would  show  that 


(Yft) 

— 

- (YH) 

.u 

n 

i.i 

n 

(vf) 



- (yh) 

n 

i j 

n 

ji 

(Yt<P) 

_ 

- (Yt<fr) 

n 

ij 

n ' 

ji 

An  efficient  method  of  computing  (31)  which  takes  (38)  and  (39)  into 
account  is  described  in  a subsequent  report. 


TV.  E-FIELD  SOLlH'f ON 

The  E-field  solution  is  obtained  by  applying  the  method  of  moments 
to  the  E-field  integral  equation.  The  E-field  integral  equation  is  de- 
rived by  setting  the  component  tangential  to  S of  the  total  electric  field 
equal  to  zero  on  S. 

The  boundary  condition  that  the  total  tangential  electric  field 
is  zero  on  S is  written  as 

- ' ES  = I K1 

ti  - tan  n -tan 
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on  S 


(AO) 


where  Rs  is  the  electric  field  due  to  the  electric  surface  current 
on  S,  R^  is  the  incident  electric,  field  given  by  either  (1)  or  (2) 
and  n is  the  intrinsic  impedance.  The  subscript  tan  denotes  tan- 
gential components  on  S.  The  1/n  terms  are  included  in  (40)  to  give 
it  the  dimensions  of  current. 

The  field  R'  can  be  expressed  in  terms  of  a vector  potential 
A and  a scalar  potential  ‘1  as 


f = - ;juA(J)  - V <t> ( J) 


where 


A(J)  = u 


„-.ik  I r-r*  I 

ds  ' 

An | r-r ' | 


*(J)  = - 
e 


JJ 

S 


e-jki'r-r'l 

Ur  | r-r ' | 


ds' 


(41) 


(42) 


(43) 


Here,  r and  r*  are  vectors  to  the  field  and  source  points  respectively 
J(r')  is  the  electric  surface  current  on  S,  k is  the  propagation  con- 
stant, u is  the  permeability,  e is  the  permittivity,  and  a is  the  sur- 
face charge  given  by 


a 


- i—  lim 

-1“’  AS  >-0 


J(r)  • u dc\ 


AS 


--  V • J(r) 


(44) 


where  u is  the  unit  vector  tangential  to  S and  normal  to  the  curve  t. 
-n 

which  bounds  the  small  portion  AS  of  S.  u points  away  f^om  AS.  The 
operator  V • is  the  surface  divergence  on  S. 

Following  the  method  of  moments,  we  write 


(I 


»,  ! 


n ] 


l' 

'll 


,h  ,( 

r .! 1 i 

n i • n 1 


I 


4- 


(45) 


wmmimwmmmm 


where  I ^ . and  1^,  are  coefficients  to  he  determined  and  J4".  and 

nj  nj  ^n.i  •*-nj 

are  expansion  functions  defined  by 


u f . (t)e 

.1 


;in<f> 


u f (t)e 
^ j 


(46) 


Next,  we  take  the  integral  over  S of  the  dot  product  of  (40)  with  each 

one  of  a collection  of  testing  functions  W*".  , . defined  by 

-mi  —mi 


Wt.  = u f.(t)e"im<l> 
— ^iru  — t l 


sL  - 


(47) 


to  obtain  the  matrix  equation 


J ([7tC]  I1  + [7^]  P)  = VC 

mn  n mn  n m 


y ([ZK]  P + [ZW]  P)  = v4, 

L mn  n inn  n m 
n 

where  the  Zfs  are  square  matrices  whose  i j-th  elements  are  defined  by 


(48) 


. > f, 

mri  i j n J J 


(7PC1). 


Wp-  * (iwA(Jq.)  + Vi(.]q.)  )ds 
'-  mi  "•  - n. i - -m 


(49) 


where  p may  he  either  t or  <t>  and  q may  he  ei  ther  t or  (J>  and  n = /p  />  is 
the  intrinsic  impedanc 
elements  are  given  by 


• > t - >7 1 f) 

the  intrinsic  impedance.  Also,  V and  V are  column  vectors  whose  i-th 
' mm 


mi 


n J 


WP  • K ds 

-mi  — 


(50) 
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where  p may  be  either  t or  <j>.  Lastly,  it4  and  are  column  vect 

. , n n 


ors 


cf  the  coefficients  T , and  1 , appearing  in  (45). 

nj  nj 

Tine  following  manipulations  ser\'e  to  transfer  the  differential 
operator  on  <J>  in  (49)  to  . Since  S is  closed. 


V * ($  W)ds  = 0 

~-S  **- 


(51) 


where  W denotes  Vr  , . Now,  the  representation 
^ -wni 


v * W = “ TZ  (P  w ‘ UJ  + 7 TT  (W  • u. ) 

•^s  ^ p 3 1 **  — t p 3(J)  ^ *ip 


(52) 


of  the  surface  divergence  and  the  definition 


„ * 3®  , 9<J> 

V $ = u — - + u — — - 

— S H ^ pd(p 


(53) 


of  the  surface  gradient  imply  that 


V • (<t>  W)  = i’V  • N + W • V i>  (54) 

<V.g  A*.  AAg 

In  (54),  the  surface  gradient  of  $ can  be  replaced  by  the  ordinary 
three-dimensional  gradient  of  <J>  because  (53)  is  the  component  of  the 
three-dimensional  gradient  tangential  to  S and  W is  tangential  to  S. 
Substitution  of  (54)  into  (51)  and  then  (51)  into  (49)  yields 

(55) 


(56) 

Since,  as  shown  in  Appendix  B,  (55)  is  zero  for  in  / n,  (48) 


(zpqi 

inn  ij 


A(.rq  ) 
- -nj 


+ 


mi 


4>(J  1 ,)  )ds 
-nj 


wiiere 


mi 


_L 

jw 


2£i 
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reduces  to 


fztc 

zt<r 

* r ‘ 

" ->t  1 

V 

n 

n 

n 

n 

z'f>t 

ZH 

-e- 

v,f> 

n 

n 

..  n j 

n 

, n=0 , +1  , -r2  , . . . 


where  is  7?”^  of  (55).  If.  is  also  shown  in  Appendix  R that  the 
n nn  1 

elements  of  Z^Q‘  are  given  bv 
n 


<zic>ij  - j 


dt.  ] dt'  [k“p f ^ (t) p ' f (t ' ) (G„sin  v sin  v'  + cos  v cos 


i;  (ofi(t)>  jf- 


(z+h., 

n ij 


dtpf^(t) 


dtf(k2p'f  ft’)  G,sin  v'  + - (p'f  .(f))C.) 

J 6 P 3t  j 4 


(Z^) 

k n ;ij 


dt  [ dt ’p  ’ f . (t ') (k2pf . (t)  0 sin  v + ~ ~~  (p f ( t) ) G ) 
J .1  1 o p dt.  i 4 


<C>ii  - j dtpfi(t) 


dt’p'f  (t')(k2n  - - — r g. ) 

.1  5 pp  A 


where  v Is  the  angle  between  the  tangent  to  the  generating  curve  and 
the  z axis  and  where 


G4  " 


17  e-1kR 


<4'  ~k¥ — COS 

„-ikR 


G,  - f d<j>  * ~~ — cos  <f ' cos  (n<*>') 
0 


- 1 kR 

G>6  = j diji 1 — — — sin  4.  * sin  (nt  ’ ) 


(57) 

v’) 

(58) 

(59) 

(60) 
(61) 

(62) 

(63) 

(6A) 


R = \/(p-k  ’ )2  + (z-z*)2  + 4po'  sin  (65) 

Hero,  p,  z , and  v depend  on  t while  p',  z',  and  v'  depend  on  t'. 

To  evaluate  (58) -(bl),  we  choose  for  p f . ( t ) the  four  impulse 
approximation  (30)  to  a triangle  function  which  reads 


pf . (t)  = — l T , <5(t-t  „ . „) 
1 k L,  p+4i-4  p+2i-2 


P=1 


(66) 


where  <S(t)  is  the  unit  impulse  function  and  T.  and  t.  are  defined  by 

(29)  and  (27)  respectively.  For  - — (pf.(t)),  we  choose  the  four  impulse 

dt  l 

approximation 


dt  (pfi(t))  I,  rp+4i-4  tp+2i-?) 


p=l 


(67) 


to  the  derivative  of  the  triangle  function  as  shown  in  Fig.  3<  Figure  4 
il 
by 


illustrates  (67).  The  coefficients  . appearing  in  (67)  arc.  given 

P ~r  +4  t — 4 


t! 


cl  2 i — 1 


4i  - 3 d„.  , + cL. 

2 i -1  > i 


tu, 

i 1 


4 i -2  iL.  , + d. 


i+1 


( 6 K) 


<i  + d 
2 l + l 2 it  2 


T 

u 1 


~l,2  i +2 

ll2  i + 1 ‘ dMt2 


.’lir  re  d . is  ile  f i Iit'il  h V ( 2 M ) . 


Sufis  I i I ul  i n i'.  ((>())  .mil  ( t)  7 ) into  ( 6 ><)  - ( U I / , wo  oh ! a i n 


2 i - I 2 i - I 1 2 i 


2 i 2i  + 1 2 i + I 2i  + 2 2i  + 2 2 i + 3 


U-rivative  of  trianple  function  (solid)  ami  four  impulse 
ipp  ro  >;  i T’l.'i  t i o;i  (arrows). 
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r )i1  = ,1  ^ ?1tTp,Tfll(r,58ln  vi'sin  V1‘  + G4cos  Vi'cos  vjt>  ■■  Tp’Tq'G4}  (69) 


(Z*c).,  - 

n l'i 


4 4 


X J^VW1"  V+5qTVW 


4 4 


(Z^).,  = y I (T  ,T  , G.sin  v , T’ ,T  ,G. } 

n 7 JO  A P q 6 i'  kp4,  p’  q 4 


P=1  q=l 


«.%-J  \ Vvvo-rr^-v 

P=1  q=l  K k p^.p  ^ t 


where 


p'  - p 4-  4i  - 4 


q'  = q + 4j  - 4 


i*  = p + 2i  - 2 


r - q + 23  - 2 


The  subscript  i'  denotes  evaluation  at  t^ , * The  subscript  j ' denotes 
evaluation  at  t ^ . In  (69)  - (72),  G^,  G,.,  ar.d  G ^ are  given  by  (62)  - (64) 
in  which  K of  (65)  is  evaluated  at  t = t^,,  t'  = t^ , which,  in  terms  of 
cylindrical  coordinates,  is  at  p,  z,p',z'  = Pj»  »Z£  * »Pj  » »Zj  » • If  i'  = .1  ’ , 
we  replace  R by  the  equivalent  distance  given  by  (36)  which,  with  i 
replaced  by  i'  reads 


(dt,/4)2  + 4p2, sin2  (>72) 


The  N point  Gaussian  quadrature  formula  (37)  which  reads 


rood*-  - f f 

k»l 


(xk  + D) 
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'■WWW 


i i.ii  1 1 


is  used  to  calculate  the  Integrals  G , G , and  G defined  by  (62)  - (64). 

4 b b 

Since  replacement  of  (i,.1,p,q)  by  (j,i,q,p)  In  (73)  implies  re- 
placement of  (1 ' , j ' ,p ' ,q  ' ) by  (j ' ,1 ' ,q ' ,p ’ ) in  (69)  - (72),  and  since 

G , G. , and  G,  are  symmetric  in  i*  and  i ' , it  is  evident:  chat 
4 b b 


(7tC) 

n ;ij 


/„tt. 

(Zn  5ji 


(zf) 


n :Lj 


- (zt4>) 

1 n 


<zn% 


(zH) 


n ji 


An  efficient  method  of  computing  (69)  - (72)  which  takes  advantage  of  (76) 
is  described  in  a subsequent  report. 


V.  COMBINED  FIELD  SOLUTION 

In  this  section,  we  assume  that  the  incident  field  (K*,H*)  due 
to  sources  outside  the  perfectly  conducting  body  whose  closed  surface  is 
S induces  a unique  electric  surface  current  J on  S and  that  this  J satis- 

A/*..  AA- 

fies  (3)  and  (40)  which  read 


- n x HS(J)  = n x H1 


just  inside  S 


- F.®  (J)  = - F1 

n -tan  - n —tan 


on  S 


s s 

where  n is  the  unit  outside  normal  vector  to  S,  (Fi  ,H  ) is  the  field  due 
to  J,  and  the  subscript  tan  denotes  tangential  components  on  S.  The  ques- 
tion arises  whether  (77)  alone  is  sufficient  to  determine  J,  whether  (78) 
alone  is  sufficient,  or  whether  information  must  lie  drawn  simultaneously 
from  both  (77)  and  (78) . 

if  both  J an J J + J satisfy  (77),  then 

H®  (J)  = 0 |ust  inside  S 

-tan  -- 


(76) 


(7  7) 
(78) 


(79) 
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Maxwell's  equations  will  he  satisfied  if 


V x y x Hs(j)  = k2HS(J)  inside  S (80) 

The  solution  J to  (77)  is  not  unique  for  values  of  k at  which  (79)  and 
(80)  admit  a nontrivial  solution  J.  This  J will  be  called  a magnetic 
cavity  mode.  Similarly  the  solution  J t:o  (78)  is  not  unique  for  values 
of  k at  which 


R'S  (J)  =0  on  S 

-tan  — 

(81) 

V x v x ES(J)  = k2ES(J) 

(82) 

admit  a nontrivial  solution  J called  an  electric  cavity  mode.  Comparison 
of  (79)  and  (80)  with  (81)  and  (82)  shows  that  the  magnetic,  and  electric 
cavity  modes  occur  at  the  same  values  of  k and  that  the  magnetic  field  of 
the  magnetic  mode  is  proportional  to  the  electric  field  of  the  electric 
mode  i-nside  8.  The  electric  mode  field  vanishes  outside  S because  its 
tangertlal  electric  field  is  zero  just  outside  S,  hut  the  magnetic  mode 
field  does  not  vanish  outside  S because  its  tangential  magnetic  field  is 
not  zero  lust  outside  S. 

The  assumed  existence  of  a unique  solution  J to  the  physical 
problem  implies  that  the  Incident  field  is  orthogonal  to  a nontrivial 
solution,  if  it  exists,  of  the  adjoint  field  problem.  The  adjoint  mag- 
netic field  operator  has  been  determined  by  Marin  [9]. 

We  will  show  that  the  solution  J to  the  combined  field  formu- 
lation [4,5] 


H (J ) - - E (J) 
- - n -tan  ~ 


x H 


~ e1 

0 -^tan 


just  inside  S 


(83) 


is  unique  and  satisfies  botli  (77)  and  (78)  whenever  a is  a positive  real 
number. 


[9]  L . Marin,  "Natural -Mode  Representation  of  Transient  Scattered  Fields," 
TREK  Trans . Antennas  Propagat . , vol . AP-2 1 , pp.  809-818,  Nov.  1973. 
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The  solution  to  (83)  is  unique  if 


- n x h8(J)  - ” E®  (J)  = 0 (84) 

^ ,vk  n *-"£  an  ~ 

implies  that  J --  0.  Scalar  multiplication  of  (84)  by  its  complex  conju- 
gate and  integration  over  S lead  ■ to 

j (ES(J)  x HS*(J))  * (-n)ds]  = 0 

(85) 


(!h°  (J) 

an 


+ -av 


E"  (J) 
1 "-tan  — - 


" ) ds  + 


2_a 

n 


[Reai. 


The  bracketed  quantity  in  (85)  is  the  real  power  fj owing  inside  S and  hence 
is  either  zero  if  the  media  is  loss-free  or  greater  than  zero  if  the  media 
is  lossy.  Thus,  (85)  implies  that 


HS  (J)  = 0 lust  inside  S (86) 

— t an  — 

ES  (J)  = 0 on  S (87) 

—tan  •— 


g 

Since  (87)  implies  that  H ^ ( J)  is  also  zero  just  outside  S,  we  obtain  the 
desired  result  that  J = 0.  Hence  the  solution  to  (83)  is  unique.  The 
statement  on  page  224  of  [5]  that  (83)  has  an  infinite  number  of  solutions 
at  eigenf requencies  is  not  correct. 

Tf  (83)  is  true,  then  (84)  - (87)  are  valid  with  HS(J)  replaced 
by  HS(J)  + and  ES(J)  replaced  by  ES(.T)  + E^.  Therefore,  (83)  implies 
noth  (77)  and  (78) . 

Since  (83)  is  the  linear  combination  of  (3)  and  (40)  with  relative 
weight  a,  the  method  of  moments  formulation  obtained  from  (83)  is  the  same 
.linear  combination  of  (17)  and  (57).  Hence, 
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where  all  matrices  and  column  vectors  have  the  same  meaning  as  in 

>(f) 

Sections  III  and  IV  except  that  I and  I now,  when  substituted  into 

' n n 

(45) , give  the  combined  field  solution  for  J. 

VI.  FAR  FIELD  MEASUREMENT  Ah'D  PLANE  WAVE  EXCITATION 

In  this  section,  measurement  vectors  are  used  to  obtain  the  far 
field  of  the  surface  current  J.  The  plane  wave  excitation  vectors  needed 
for  the  H-field,  E-field,  and  combined  field  solutions  are  then  expressed 
in  terms  of  these  measurement  vectors. 

By  reciprocity. 


(J) 


K(Ii)ds 


(89) 


where  ES(J)  is  the  far  electric  field  due  to  J,  I?,  is  a 

dipole  at  the  far  field  measurement  point,  and  E(I£^)  is 

due  to  I£  evaluated  at  r on  S . If  £ is  tangent  to  the 
— r — r 


receiving  electric 
the  electric  field 
radiation  sphere, 


E(i*r) 
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-jkr 


4ur 


I £ 
-r 
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(90) 


where  r is  the  distance  between  the  measurement  point  and  the  origin  in 
the  vicinity  of  S and  is  the  propagation  vector  of  the  plane  wave  coming 
frc-m  I£_.  Subs t i ting  (45),  (46),  and  (90)  into  (89)  and  letting  £^  be 

IT  f 

either  u"  or  u.,  we  obtain 
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given  by 
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dtp  (t) 
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u1)e 


K-k 


r + n(<j>-<t>r)) 


(92) 


where  <j>  is  the  azimuth  of  the  far  field  measurement  point.  In  (91),  Eg(J) 
and  F,S(J)  are  respectively  the  and  u*'  components  of  ES(J). 

(p  ■'wtp  <w  *** 

With  a view  toward  evaluation  of  (92),  we  note  from  Figs.  2,  A-l, 
and  A~2  that 


* • >4  ■ 

- sin  6 cos  v + cos  0 sin  v cos(d)-<t>  ) 
r r r 

’ 'if) 

- cos  0 sin (<f>-i>  ) 

r r 

(93) 

r 

~t  ' ^ 

sin  v sin(<t>~<}>r) 

~<f>  ' ^ 

cos  (4>—<fiT.) 

-k  • r = 

~T  - 

kz  cos  0 + kp  sin  9 cos(d>-d>  ) 

r r r 

Substituting  (93)  and  (30)  into  (92)  and  taking  advantage  of  the  integral 
formula 

2 TT 

-n  r j (kp  sin  0 cos  cf>  + n<j>) 

J (kp  sin  0 ) = 4; — I e r dd>  (94) 

n r l w j 

0 

deduced  from  (9.1.21)  of  [10]  for  Bessel  functions,  we  obtain 


[10]  M.  Abramowitz  and  T,  A.  Stegun,  "Handbook  of  Mathematical  Functions," 
U.  S.  Government  Printing  Office,  Washington,  D.C.  (Natl.  Bur.  Std. 

11.  S.  App  1 . Math.  Scr.  55),  1964,  p.  360. 
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(95) 
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where 


J ~ J (kp  sin  0 ) 
n n r 


(96) 


In  (95),  p,  z y and  v are  to  be  evaluated  at  t = t 

p+2i~2  * 

Substitute  (1)  and  (2)  into  (18)  and  (19)  to  obtain 
2tt 


fpq  = 

ni 


dtpfjCt) 
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t -1  * I - n^) 

d cb ( u x n)  • (k  x u ) e 
~P  " ~t  ~q 


(97) 


and  then  substitute  (1)  and  (2)  into  (50)  where  WP  is  given  lvr  (47)  to 
ob  tain 


pq  f f t j(--k.  • r - n<f>) 

^ -j  ~ k | dtpf  (t)  d <(>  ( u • u )e 

n 1 J 1 j ~p  ~q 


(98) 


where  p is  either  t or  <f>.  Tlie  additional  superscript  q on  (lie  left-hand 
sides  of  (97)  and  (98)  is  either  0 or  <J>  according  as  the  incident:  electric 
field  is  0 polarized  as  in  (1)  or  4,  polarized  as  in  (2).  Comparison  of 
(97)  with  (92)  shows  that 
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where  the  R's  on  the  right-hand  side  of  (99)  are  to  be  evaluated  at 
er  = • Comparison  of  (98)  with  (92)  leads  to 

VPq  = RPq 
n -n 


(100) 


where  p is  either  t or  <p,  q is  either  6 or  <J>,  and  Rpq  is  to  be  evaluated 

-n 

at  0 = 0 , 

r t 

The  expressions  (45)  and  (91)  for  the  electric  surface  current  and 
far  field  can  be  simplified  by  combining  the  +n  and  ~n  terms.  Substitu- 
ting (46)  into  (45) , we  obtain 


f . I {(f  T»)u  + (f  !*<),.  > (101) 

ns3 oo  n n 

where  f is  the  transpose  of  the  column  vecLor  f of  the  f.(t),  and  ?tq  and 

I^q  are  column  vectors  of  the  coefficients  and  1^ , respectively.  The 
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additional  superscript  q is  either  0 or  if  according  as  the  incident  elec- 
tric field  is  6 polarized  as  in  (1)  or  <J>  polarized  as  in  (2).  The  column 
-*tq  -y-rhn 

vectors  1^  and  I appearing  in  (101)  are  obtained  by  solving  either  the 
H-field  matrix  equation  (17),  the  E-ficld  matrix  equation  (57),  or  the 
combined  field  matrix  equation  (88)  with  the  additional  superscript  q on 
the  column  vectors  therein  to  denote  the  polarization  of  the  incident 
electric  field. 

Inspection  of  (20)  - (26)  and  (58)  - (64)  reveals  that 
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and 


» n=0»l,2. 


From  (95),  it  is  apparent  that 
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Substitution  of  (104)  into  (99)  and  (100)  gives 
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Since  the  properties  (102)  and  (103)  survive  matrix  inversion,  it  is  evi- 
dent from  (17),  (57),  (88),  (105),  and  (106)  that 
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Fn  view  of  (107),  (101)  becomes 
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= (f  1^)  u + Y (21(f  Ttt|>)u  sin  (nip)  + 2(f  I^)u  cos(n<f)} 
~ o ~<b  u , • n ~t  n -~<j> 

n=  I 


where  f is  a row  vector  of  the  (t)  . If  pf^(t)  is  the  triangle  function 
itself  rather  than  the  four  impulse  approximation  (30)  to  the  triangle 
function,  then 
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Equations  (104.1  and  (107)  reduce  the  specializations  of  (91)  for 
0 and  <f>  polarized  Incident  electric  fields  to 
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In  (110),  the  first  subscript  on  E denotes  the  component  under  con- 
sideration whereas  the  second  subscript  denotes  the  polarization  of 
the  incident  electric  field. 


The  scattering  cross  section  o is  that  area  for  which  the  incident 
wave  contains  enough  power  to  produce  by  omnidirectional  radiation  the 
scattered  power  density  at  the  far  field  point.  Specializing  a to  the 
four  different  polarizations,  we  obtain 


4tr' 


2nc  F.  (J) 



(111) 


where  p is  either  0 or  and  q is  either  6 or  4) . In  (111),  p is  the 
receiver  polarization,  q is  the  transmitter  polarization,  X = 2ir/k  is 

g 

the  wavelength,  and  E (J)  is  given  by  (110). 
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VI T ■ EXAMPLES 

Computer  program  subroutines  have  been  written  to  calculate  the 
square  matrices  and  measurement  vectors  needed  for  the  K-field,  E-Field 
and  combined  field  solutions.  These  subroutines  will  be  described  and 
listed  in  a subsequent  report.  Some  computational  results  obtained  with 
these  subroutines  are  given  in  this  section. 


We  compare  computed  approximations  to  the  electric  current  induced 
on  tlie  surface  of  a conducting  sphere  by  an  axial  incident  plane  wave  with 
the  known  exact  solution.  Figure  5 is  a plot  of  A versus  ka  where  k is 
the  propagation  constant,  a is  the  radius  of  the  sphere,  and  A is  defined 
by 
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V lwF7/T»_ 

id 


(112) 


Here,  S is  the  surface  of  the  sphere  and  11*  is  the  'incident  magnet  i c 
field.  Also,  V is  the  exact  electric  cut  rent  given  by  (6-101)  of  [61 
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Fleet rir  nirmit  error 
excited  by  an  i nei  dent 
r i roles  denote  K-t  ield 
solution.  Values  ol  A 


A versus  ka  for  a eondue  t i n p,  sphere  of  radius 
plant*  wave.  Squares  denote  11-field  solution, 
solution,  and  triangles  denote  combined  field 
■ 0.5  are  plot  t ed  at  t lie  t op  . 


a 


q 

and  J is  the  computed  approximation  to  J . The  squares,  circles,  and 
triangles  in  Fig.  5 tell  which  field  solution  (E-field,  H-field,  or 
combined  field  with  a = 1 in  (88))  J is  obtained  from.  In  (112),  the 
expression 

|j  - JC|2  = (J  - Jc)  • (J  - jf)*  (113) 

where  * denotes  complex  conjugate  is  the  time  average  of  the  square  of 
the  length  of  the  time  dependent  vector  whose  root -me an -square  phasor 
is  (J-Jc) . 

The  number  P (maximum  value  of  (i+1)  in  (27))  of  data  points  on 
the  generating  curve  is  31  for  all  the  examples  of  this  section.  These 
data  points,  equally  spaced  from  the  lower  pole  to  the  upper  pole  of  the 
sphere,  give  14  expansion  functions  for  the  t directed  electric  current 
and  14  expansion  functions  for  the  <f>  directed  current.  Both  surface 
integrals  in  (112)  are  evaluated  by  integrating  analytically  in  <j>  and 
by  sampling  in  t at  the  30  points  defined  by  (27).  The  number  N of 
points  used  in  the  Gaussian  quadrature  integration  (37)  is  20. 

Tn  Fig.  5,  the  H-field  solution  is  generally  the  best  and  the 
F-field  solution  is  the  worst  as  far  as  J is  concerned.  For  both  F.-field 
and  H-field  solutions,  the  error  in  .J  is  large  at  resonances  of  the  spheri- 
cal cavity  which  are  tabulated  on  page  270  of  [6].  The  combined  field  solu- 
tion is  not  affected  by  these  resonances.  Note  that  the  error  in  J for  the 
F.-fleld  solution  has  a peak  around  ka  = 1.35  which  is  far  from  any  resonance 
This  peak  disappeared  when  the  number  P of  data  points  was  reduced  from  31 
to  21. 

, 2 

Figure  6 compares  the  noimalized  radar  cross  section  o/ita  in  the 

backscattering  direction  obtained  from  the  H-field  solution,  the  K-fielc 

solution,  and  the  combined  field  solution  with  the  exact  a/ The 
2 

exact  o/ira*"  is  calculated  from  (6-105)  of  [6]  in  which  A is  our  o.  In 

e 

Fig.  6,  the  squares  denote  the  H-field  solution,  the  circles  denote  t he 
E-field  solution,  the  triangles  denote  the  combined  field  solution,  and 
the  solid  line  denotes  the  exact  solution.  The  E— field,  H-field  and 
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combined  field  solutions  for  a flux'  arc  beginning  to  deteriorate  for  the 
larger  val  ies  of  ka  in  Fig  6.  At  ka  = 6 » there  are  only  5 expansion 
function.-;  per  wavelength. 

Figures  7 and  8 show  the  curves  of  Figs.  5 and  6 respectively 
in  more  detail  in  the  vicinity  of  the  first  resonance  which  occurs  at 
ka  - 2.744.  The  disturbance  in  the  11-field  solution  occurs  quite  close 
to  the  resonant  frequency,  hut  the  disturbance  in  the  E-field  solution 
occurs  at  a slightly  higher  frequency.  Although  the  error  in  .1  for  the 

■A 

F.-f leld  solution  is  tremendous,  the  error  in  c/tux'  for  the  F.-field  solu- 
tion is  quite  small  except  at  two  or  three  points.  The  following  expla- 
nation is  offered.  According  to  Section  V,  the  F- fie Id  solution  does  not 
determine  how  much  of  the  electric  cavity  mode  current  is  contained  in  J. 
Hence,  we  suspect  that  our  numerical  F- field  solution  does  not  contain 
the  right  amount  of  the  electric  cavity  mode.  If  this  suspicion  is  true, 
then  the  radar  cross  section  can  still,  he  quite  accurate  because  the 
electric  cavity  mode  does  not  radiate  any  external  field. 

Figures  7 and  8 show  that  the  combined  field  solution  is  much 
better  than  either  the  h— field  solution  or  the  K~fie1d  solution  in  the 
vicinity  of  the  first  resonance. 
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An  H-fie'id  solution,  an  E-field  solution  and  a combined  field 
solution  for  plane  wave  scattering  from  a perfectly  conducting  body  of 
revolution  have  beet)  developed  and  expounded.  The  H-f leld  solution  is 
a mod  i f f cat  ion  and  general  izat  i >n  to  oblique  incidence  of  'ds ! engh  i ' s j 2 j 
i tnp ..  i »e  solution.  Instead  of  impulse  expansion  fine*  ions  and  f imp  son  ' s 
rule1  tor  i n fog  rat  i on  with  respect  to  the  azimuth  : , we  use  four  impulse 
npprc»:  i mat  i ous  to  triangle  I unctions  and  Causs  tan  quaelrate.ro  integration. 
The  ''.-field  solution  is  t hat  ol  [ 1 ] with  impulse  (1  roan's  f unct  ions  obtained 
I com  Eaussian  quadrature1  integration  in  the  azimuth  ■!  instead  e > f pulse 
Ere  n's  fun<-t  ions  oh.  aimsl  t rmn  equal  interval  , equal  weight  samp  ling.  Tin1 
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ig.  7.  Fleet  rie  current  error  A versus  ka  for  a conducting  s 
a excited  !>v  an  incident,  plane  wave.  Squares  denote 
circles  denote  K- field  solution,  and  triamjes  denote 
solution.  The  first  resonjii.ee  of  the  splierieal  eavi  t 
Values  of  A > 0.5  are  plotted  at  the  top. 


sphere  of  radius 
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Fig.  8.  Normalized  radar  cross  section  o/n a'  versus  ka  for  a conducting  | 

sphere  of  radius  a.  Solid  line  denotes  exact,  solution,  squares  ] 

denote  li-field  solution,  circles  denote  K-fleld  solut  ion,  and  j 

triangles  denote  combined  field  solution.  Tin*  first  resonance  1 

of  t!u>  spherical  cavity  is  at  ka  ~ 2 . 7 44 . 'j 


combined  field  formulation  is  a linear  combination  of  t lie  conations  for  the 
H- field  and  F.-f ield  solutions.  Computer  program  subroutines  which  calculate 
tlie  square  matrices  for  the  H-field  and  I'-field  solutions  and  the  plane  wave 
measurement  vectors  will  appear  in  a f ortheoming  report  . Tlie  subroutine  which 
calculates  the  square  matrix  for  the  '.-field  solution  executes  appreciably 
faster  and  requires  considerably  less  storage  to  obtain  the  same  kind  of 
accuracy  as  a previous  solution  (ll|. 

The  H-field  and  F.-f  it:  Id  solutions  deteriorate  in  ‘lie  vicinity  of 
cavity  resonances  because  their  homogeneous  equations  admit  nontrivial  solu- 
tions at  these  resonances.  In  Section  V,  it;  is  shown  that  t he  homogeneous 
equation  associated  with  the.  combined  field  formulation  has  no  nontrivial 
solution  if  the  relative  weight  n of  the  F.-f  ield  equation  is  real  and,  if 
the  inside  media  is  lossy,  positive.  For  this  reason  the  combined  field 
solution  is  much  better  than  either  the  H-field  or  R-field  solutions  in  the 
vicinity  of  cavity  resonances.  Figures  7 and  8 bear  this  out 

For  the  examples  of  Section  VII,  the  relative  weight  u of  the  F, -field 
equation  in  the  combined  field  formulation  is  unity.  This  puts  the  H-field 
and  F.-f ielcl  equations  on  a more  or  less  equal  footing  because  the  magnitude 
of  the  excitation  due  to  the  H-field  equation  is  then  that  due  to  the  E-field 
equation  rotated  90°  in  space.  According  to  vig.  6,  the  H-field  solution 
for  ct/tt a is  generally  a bit  more  accurate  than  the  F-fielrJ  solution  away 
from  the  cavltv  resonances.  This  suggests  that  one  weights  the  K-field 
equation  less  than  the  H-field  equation  in  the  combined  field  formulation. 

Oshiro  et  a.l . [4,5]  conclude  from  their  plots  of  mean  error  versus  a 
for  0 < rx  <_  1 that  an  <*  value  on  the  order  of  0.2  is  best.  However,  there 
is  no  logical  reason  for  ruling  out  negative  values  of  a when  the  electric 
surface  current  radiates  into  a loss-free  inside  media  because  then  the  left- 
hand  side  of  (85)  does  not  depend  on  the  sign  of  a.  We  see  little  significance 
in  the  facts  that  the  magnitude  of  the  combined  field  excitation  on  the  right- 
hand  side  of  (83)  is  generally  larger  on  Lire  illuminated  poition  of  the  surface 
of  the  body  of  revolution  than  in  the  shadow  zone  for  a 0 and  that  the  oppo- 
site is  true  for  <t  - 0. 

I ’•  1 ] F . !•  . Harrington  and  .1 . R.  Maurz,  "Radiation  and  Scattering  from  Bodies 
of  Revolution,"  Report  AFCK1, -69-0 305 , Contract  No.  F-] 9628-67-0-0233 
between  Syracuse  University  and  Air  Force  Cambridge  Research  laboratories, 
'ul v 1969. 


APPENDIX  A 


DERIVATION  OF  THE  H -FIELD  INTEGRAL  EQUATION 


The  purpose  of  Appendix  A is  to  obtain  (9)  from  (7)  and  (8) . In 
view  of  (8) , 


(r-r')  x j(r')  - (r-r')  x u’  Jt(t\<f>')  + (r-r')  x u’  JT  (t ' , $ ' ) (A-l) 


The  cross  products  on  the  right-hand  side  of  (A-l)  are  evaluated  by  expres- 
sing all  vectors  in  term; 
directions  respectively. 


sing  all  vectors  in  terms  of  unit  vectors  u , and  u in  the  p,  and  z 

~p  -(f)  ~~z 


r = 

up  +•  u z 
~P  -z 

(A-2) 

r ' = 

U^p  'cOs(4>  ' -<)>)  + 

u p ' si n(4> '--())) 

+ u z ' 

— z 

(A-3) 

Hj- 

u sin  v'  cos(<j>'- 

~p 

-cp)  + u sin  v' 

sin(<)> ' -4O  + u^cos  v' 

(A-4) 

yj= 

- u sin  (A  ' — 4> ) + 

-p 

U , COS  (<P  '-<)>) 

(A-5) 

Equation  (A-3)  has  been  obtained  bv  first  writing 

r ' = u'o  ' + u ?.'  (A-6) 

— —p'  — z 

and  then  using  Fig.  A-l  to  express  u’  in  terms  of  u and  ui . To  verifv 

~P  -p  -N> 

(A-4) , use  Fig.  A-2  to  express  u*  in  terms  of  u'  and  u and  then  use 

~ t -p  —z 

Fig.  A-l  to  express  u'  in  terms  of  u and  u., 

~P  ~P  ~<t> 

Substitution  of  (A-/)  - (A-5)  into  (A-l)  yields 

(r-r')  x J(i')  - (u  (-p'cos  v'  + (z'-z)  sin  v')  sin(<h'-d>) 

- “ ~ “ -P 

+ n (-(p-p 'cos (<(>'-<}>))  cos  v'  - (z'-z)  sin  v'  cos (<j> ' -tp) ) 

+ u p sin  v ' sin  (<J>  ' — <h)  } It(t  ’,<h')  + (u  (z'-z)  cos  (<{>  * -<j>) 

-z  ~P 

+ u (z'-z)  sin  (<)' ' -<j>)  i-  c ' + p cos  ( ' --<!)  ) } .T^(t ' 1 A - 7 ) 

" £ 
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To  find  the  u and  u components  of  (7) , we  need 

*~t  ■"•p 


it  * n x [(r-r')  x J(r')]  ~ u • f(r-r')  x J(r*)] 


u • n x [ (r-r * ) x J(r' ) ] = - u • [ (r-r ' ) x J(r ') ] 


With  the  help  of 


u = u sin  v + u cos  v 
— t — p ~z 


and  (A-7) , (A-8)  and  (A-9)  become 


(A-8) 

(A-9) 


(A-10) 


u • n x [(r-r')  x J(r_' )]  = {((p'-p)cos  v'  - (z'-z)  sin  v' ) cos ($ ' -9) 


- 2p  cos  v'sin2(^~^)}  J^t’,*')  + (z'-z)  sin(9  •-$)  J*(t  * , <fi  ’ ) 


(A-ll) 


• n x [(r-r')  x J (r ' ) ] = (p'sir.  v cos  v'  - p sin  v'  cos  v 

- (z'-z)  sin  v sin  v')  sin((j>  '-<(>)  .Jt  (t  ’ , <f> ' ) + {((p'-p)cos  v 

- (z'-z)  sin  v)  cos(<fi'-9)  + 2p’cos  v sin2^-^)}  J^(t ',<(>' ) (A-12) 

The  distance  |r-r'|  appearing  in  (7)  is  the  square  root  of  the 
sum  of  the  squares  of  (z-z’)  and  the  projection  of  (r-rj)  in  the  xy 
plane.  Hence, 


r-r 


y(z-z’) 


+ P ' 2+p2-2pp  ' cos  (9 '-<}>)  ="[/(p-p  ' ) + (z-z  ' ) 2H-4pp  ' sin" 


(A-13) 


/m  integral  with  respect  to  9 ' results  when  the  surface  integral  in 
(7)  is  Iterated.  Because  this  integral  with  respect  to  <[>'  is  an  integral 
of  a 2 tt  periodic  function  of  <J>'  over  the  period  2tr,  9'  may  be  replaced  by 
<f>*  + <f>  without  changing  the  value  of  the  integral.  Substitution  of  (A-ll)  - 
(A-13)  into  (7)  leads  to  the  desired  H-field  integral  equation  (9). 
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APPENDIX  B 

DERIVATION  OF  THE  E-FIELD  MATRIX  EQUATION 

The  testing  functions  wj\  appearing  in  (55)  are  defined  by  (47) 
From  (55),  (47),  and  (52),  we  obtain 


'mi 


A = 

mi 


JL  f (t)e->* 

ap  i 


(B-l) 


The  vector  and  scalar  potentials  A and  $ appearing  in  (55)  are  given  by 

(42)  and  (43)  with  defined  by  (46).  In  agreement  with  (44),  (46), 
and  (52),  and  in  analogy  with  (B-l),  the  charge  density  a appearing  in 

(43)  is  specialized  to  either  a1" . or  cA  given  by 

no  nj  7 


Vi  = AA  dh  (P,fj(t'))ein^ 


(B-2) 


-n 


^ f.(t’)ejnc 


n.l  mp  j 

In  view  of  the  above  considerations,  (55)  can  be  rewritten  as 


(ztc) 

mn  i ) 


(ZK)  . . 

mn  ij 


(A)., 

mn  i j 


mn  1 j 


dt 


dt’ 


2r  2rTr  -ik]  r-r ' ! j (n<i’ -m<()) 

» 6 


d<I 

0 0 


d4> ' 


^TFFT“ 


(cont.  on  next  page) 
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jkpfi(t)p'fj(t,)ut  * (Pf.(t))  gfr  (P'fjtt’)) 


jkpf1(t)p'f.(t-)u^  (p'fjCt1)) 

jkpfi(t)p,f.(t')ut  • +~r^t-  (pfi(t))p’fi(t’) 


(B-3) 


iron. 


ipf . (t)p 'f . ( t ’ ) (ku  • u' 

x n '“4>  n)>  kpp 


r) 


To  facilitate  evaluation  of  the  dot  products  appearing  in  (B-3),  we 


write 


u = u sin  v + u cos  v 
~t  "p  •"■z 


u'  = u'  sin  v'  + u cos  v' 
-t  ~p  ~z 


(B-4) 


where  u , u',  and  v'  are  defined  in  Figs.  A-l  and  A-2 . With  the  help  of 

•PD  -"p 

(B-4) , we  obtain 


u • u'  = sin  v sin  v'  cos(d>'-d>)  + cos  v cos  v' 

-t  Y 


u • u'  = sin  v'  sin(4> ' -if) 

~<j>  -t 

(B-5) 

u • u'  = - sin  v sin  (if' -<}>) 

--<p 


u , • u'  = cos  (<j)  *-<)>) 


The  distance  |r-r'|  is  given  by  (A-13)  which  reads 


| r-r' | = \j  (p-p  ' ) 2 + (z-z’)i  + 4pp'sin2  (-—^ 


(B-6) 


In  view  of  (R-5)  and  (B-6) , the  integrands  of  (B-3)  are  periodic 
functions  of  <f>  ’ with  period  2tt.  Hencey  <p ' can  be  replaced  by  (<t>'-f<J)) 
without  changing  the  values  of  any  of  the  integrals.  When  this  is  done, 
the  4>  dependence  of  the  integrands  becomes  e^n  which  when  integrated 
gives  2 7 r for  m=n  and  zero  for  m 4 n.  Taking  the  liberty  of  replacing  the 
double  subscript  mn  by  the  single  subscript  n,  we  obtain  (53)  - (61).  The 
forms  (62)  - (64)  of  the  <b'  integrals  follow  from  the  even  or  odd  symmetry 
of  the  terms  in  (B-5)  about  (4>'-<j>)  = 0 and  the  even  symmetry  of  (B-6)  about 
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